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' We discuss a sum rule satisfied by the correlation function of two particles with small relative 

momenta. The sum rule, which results from the completeness condition of the quantum states of the 
two particles, is derived and checked to see how it works in practice. The sum rule is shown to be 
trivially satisfied by free particle pairs, and then three different systems of interacting particles are 

, considered. We discuss neutron and proton pairs in the s— wave approximation and the case of the 

so-called hard spheres with the phase shifts taken into account up to I — 4. Finally, the Coulomb 
pL^ ' system of two charged particles is analyzed. 
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\Q , I. INTRODUCTION 

O 



The correlation functions of two identical or non-identical particles with 'small' relative momenta have been ex- 
tensively studied in nuclear collisions for bombarding energies from tens of MeV [1] to hundreds of GeV [2]. These 
f^S 1 functions provide information about space-time characteristics of particle sources in the collisions. As shown by one 
of us [3] , the correlation function integrated over particle relative momentum satisfies a simple and exact relation due 
to the completeness of the particle quantum states. The relation can be used to get a particle phase-space density, 
■ following the method [4,5], with no need to extract the Coulomb interaction. The sum-rule offers a constraint for 
the procedure of imaging [6,7] which inverts the correlation function and provides the source function. The sum-rule 
is also helpful to handle the correlation functions of exotic systems like pA [8] when the interparticle correlation is 
poorly known. 

The aim of this paper is to discuss the sum rule in detail, to prove or disprove its usefulness in the experimental 
studies. Therefore, we first derive the sum rule and show that it is trivially satisfied by free particles. Then, we 
analyze the correlations in the neutron-proton (n-p) system where there are both the two-particle scattering states 
and a bound state i.e. a deuteron. We prove that in spite of the attractive interaction the n-p correlation can be 
negative due to the deuteron formation. Although some qualitative features dictated by the sum rule are certainly 
seen, the calculated correlation function does not satisfy the sum rule. This is not surprising as not only the momenta 
of order 1/R, where R is the source size, but also larger momenta contribute to the sum-rule integral. Consequently, 
the s— wave approximation, which is used compute the n-p correlation function, is insufficient. 

To see how important are the higher I phase shifts, we discuss a system of the so-called hard spheres and compute 
the correlation function of identical and non-identical particles, taking into account the phase shifts up to I — 4. 
Paradoxically, the higher I contributions do not improve the situation but make it even worse. 

Finally, we discuss the correlations due to the Coulomb interaction. This is the case of particular interest as one 
usually measures the correlations functions of charge particles which experience the Coulomb interaction. And the 
integral, which is controlled by the sum rule, is used to determine the particle phase-space density [4,5]. In the case 
of Coulomb interactions, the exact wave functions are known, and consequently, the exact correlation functions can 
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be computed. Unfortunately, as we discuss in detail, the integral of interest appears to be divergent, and the sum 
rule does not hold for this most important case. 



II. PRELIMINARIES 

To avoid unnecessary complications, our considerations are non-relativistic and we start with the formula repeatedly 
discussed in the literature which expresses the correlation function R(q) of two particles with the relative momentum 
q as 

R(cO = J d 3 rV r (r)\^(r)\ 2 , (1) 

where </> q (r) is the wave function of relative motion of the two particles and T> r (r) is the effective source function 
defined through the probability density D r (r,t) to emit the two particles at the relative distance r and the time 
difference t as 

D r (r) = J dtD r (r-vt,t) , (2) 

with v being the particle pair velocity with respect to the source; the relative source function D r (r,t) is given by 
the single-particle source function D(r,t), which describes the space-time emission points of a single particle, in the 
following way 

D r (r, t)= J d 3 RdT D(R — ^r, T — ±t) D(R + ^r, T + ±t) . 

We note that T> r (r), D r (r,t) and D(r, t) are all normalized as 

J d 3 r V r (r) = J d 3 r dt D r (r, t) = J d 3 r dt D r (r, t) = 1 . (3) 

We also observe that the spherically symmetric single-particle source function D(r,t) provides, in general, the 
effective T> r (r) which is elongated along the velocity v. To simplify calculations, we, however, assume here that the 
source function T> r (r) is spherically symmetric. Such an assumption makes sense when the single-particle source 
function D(r,t) is spherically symmetric and particles are emitted instantaneously i.e. D(r,t) = D(\r\)S(t — to)- 
Then, 

V r (r) = J d 3 R D(R - ir) £>(R + ^r) . 
The single-particle source function is often chosen in the Gaussian form 

It gives the mean radius squared of a source equal (r 2 ) = 3rg, and it leads to the effective relative source function as 

^ (r) ^(4^7l cxp (-^)- ^ 

III. SUM RULE 

Let us consider the correlation function integrated over the relative momentum. Since i?(q) — > 1 when q — ► oo, we 
rather discuss the integral of i?(q) — 1. Using Eq. (1) and taking into account the normalization condition of T> r (r) 
(3), one finds, after changing order of the r— and q— integration, the expression 

/ §f (* (q) - 1 >S d3r Vrir) I ( |0q(r)|2 - x ) • (5) 
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It appears that the integral over q in the r.h.s. of Eq. (5) is determined by the quantum-mechanical completeness 
condition. Indeed, the wave functions satisfy the well-known closure relation 

/ (0 + E MWW) = ^ - r ') ± <5(3) ( r + r ') ' 

where </> Q represents a possible bound state of the two particles of interest. When the particles are not identical the 
second term in the r.h.s. of Eq. (6) is not present. This term guarantees the right symmetry for both sides of the 
equation for the case of identical particles. The upper sign is for bosons while the lower one for fermions. The wave 
function of identical bosons (fermions) </> q (r) is (anti-)symmetric when r — > r, and the r.h.s of Eq. (6) is indeed 
(anti-)symmetric when r — > — r or r' — > — r'. If the particles of interest carry spin, the summation over spin degrees 
of freedom in the l.h.s. of Eq. (6) is implied. When the integral representation of 5^(r — r') is used, the relation (5) 
can be rewritten as 



/ 



d 3 c 



(r)^(r')- e lq(r - r,) ) +^0 Q (r)C(r')=±^ (3) (r + r') . 



(2tt)3 

Now, we take the limit r' — > r and get the relation 

/ (0 ('^ (r)|2 = ± E \Mr)\ 2 ■ (7) 

When Eq. (7) is substituted into Eq. (5), one finds the desired sum rule 

f d 3 q (i?(q) - l) = ±tt 3 V r (0) - A <* ■ (8) 
where A a is the formation rate of a bound state a 

A a = (2nf J d 3 rV r (r)\Mr)\ 2 ■ (9) 

A a relates the cross section to produce the bound state a with the momentum P to the cross section to produce the 
two particles with the momenta P/2 as 

da a da 

= A n 



dP ° d(P/2)d(P/2) ' 

The tilde means that the short range correlations are removed from the two-particle cross section which is usually 
taken as a product of the single-particle cross sections. 

The completeness condition is, obviously, valid for any inter-particle interaction. It is also valid when the pair of 
particles interact with the time-independent external field, e.g. the Coulomb field generated by the particle source. 
Thus, the sum rule (8) holds under very general conditions as long as the basic formula (1) is justified, in particular 
as long as the source function T> r (r) is q— independent and spin independent. The validity of these assumptions can 
be only tested within a microscopic model of nucleus- nucleus collision which properly describes the quantum particle 
correlations and bound state formation. 

There are several potential applications of the relation (8). As already mentioned, the integral of the correlation 
function in the l.h.s. of Eq. (8) is used to determine a particle phase-space density. The method [4,5] assumes that the 
correlation function represents non-interacting particles. Consequently, before computing the integral, the correlation 
function is corrected in such a way that the Coulomb interaction is removed. However, the sum rule (8) shows that 
the integral of the correlation function is independent of the inter-particle interaction. Therefore, there is no need to 
extract the Coulomb interaction. We note that this procedure is rather model dependent. 

To get information about the source function D(r,t), one usually parameterizes the function, computes the corre- 
lation function and compares it with experimental data. The method of imaging [6,7] provides the source function 
directly, inverting the functional i?[2? r ] (1) with the experimental correlation function as an input. As seen, the 
relation (8), which gives T> r (0), can be treated as a useful constraint of the imaging method. 

The sum-rule is also helpful to better understand the correlation functions. In particular, the sum rule shows that 
the correlation function can be negative in spite of an attractive inter-particle interaction. This happens when the 
particles form bound states represented by the second term in the r.h.s. of Eq. (8) or the interaction is strongly 
absorptive as in the case of antiproton-lambda system which annihilates into mesons. Then, the scattering states of 
the pA do not give a complete set of quantum states of the system and the integral from the l.h.s. of Eq. (8) has to 
be negative. The recently measured correlation function of p and A [8] is indeed negative at small relative momenta. 
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IV. FREE PARTICLES 



In the case of non-interacting particles, the correlation function differs from unity only for identical particles. Then, 
the wave function, which enters the correlation function (1), is an (anti-)symmetrized plane wave 

= 71 ( e<qp±e-<qp ) < 

with the upper sign for bosons and lower for fermions. Then, the integration over q in Eq. (5) can be explicitly 
performed without reference to the completeness condition (7), and one finds 

fd 3 q(R(< l )-l)=±7r 3 V r (0). (10) 

It this way the sum rule (10) was found in [9], see also [4]. 

Although the sum rule (8) assumes the integration up to the infinite momentum, one expects that the integral in 
Eq. (8) saturates at sufficiently large q. To discuss the problem quantitatively, we define the function 

%max) - 4^ /"*"* dq q 2 (R(q) - 1) . (11) 
Jo 

As already mentioned, the source function is assumed to be spherically symmetric, and consequently the correlation 
function depends only on q = |q|. Thus, the angular integration is trivially performed. 

For the Gaussian effective source (4), when ir 3 V r (0) = ( v / 7r/2ro) 3 , the free functions R(q) and S(q max ) equal 

R(q) = 1 ± e- 4r « q2 , 



S{q max )=±(^j E 2/3 ((2r g max ) 3 ) 



where 



and T(z) is the Euler gamma function. Since for large x, the function E„(x) can be expressed as, see e.g. [10], 

E ^^-f(l^^( 1 + (^))' 

we have the approximation 

S(_)- ± (^)'(l-*^=^i.), 

for (2r (/max) 3 ^> 1- As seen, the sum-rule integral (10) is saturated for g max not much exceeding l/2r , and obviously, 
the sum rule is satisfied. 



V. NEUTRON-PROTON SYSTEM 



In this section, we discuss the interacting neutron-proton pair which is either in the spin singlet or triplet state. 
The nucleons, produced in high-energy nuclear collisions, are usually assumed to be unpolarized, and one considers 
the spin-averaged correlation function R which is a sum of the singlet and triplet correlation functions R s,t with 
the weight coefficients 1/4 and 3/4, respectively. Here, we consider, however, the singlet and the triplet correlation 
functions separately. Then, the sum rule (8) reads 
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j d 3 q (i? s (q) - l) = , (12) 

J d 3 q (i?*(q) - l) = -A d . (13) 

Following [11], we calculate the correlation functions R s,t assuming that the source radius is significantly larger 
than the n-p interaction range. Then, the wave function of the n-p pair (in a scattering state) can be approximated 
by its asymptotic form 

iqr 

^(r)=e""- + /*'*(<?)— , (14) 
where / s, *(q) is the scattering amplitude. It is chosen as 

^ ^ 1 — i<i s '*a s '*g 2 + iqa s ^ ' ^ ^ 

where a s,t (d 6 '*) is the scattering length (effective range) of the n-p scattering; a s = —23.7 fm, d s = 2.7 fm and 
a* = 5.4 fm, rf* = 1.7 fm [12]. The amplitude (15) takes into account only the s— wave scattering. This is justified as 
long as only 'small' relative momenta are considered. 

Substituting the wave function (14) into the formula (1) with the source function (4), we get 

RS'Hq) = 1 + Re[/'-*(g)] -L e~^W CT&{ 2r q) Im [/'•*(«)] -L(l - e~ ir W) + |/'-*( 9 )| 2 ^ , (16) 

zr o1 zr o" z 'o 



where 



2 f x 2 
erfi(a;) = — / die* 2 
V 71 " Jo 



Since the source described by the formula (4) is spherically symmetric, the correlation function (16) does not depend 
on q but on q only. 

In Figs. 1 and 2 we show, respectively, the singlet and triplet correlation functions computed for three values of 
the source size parameter r . As seen, the triplet correlation is negative in spite of the attractive neutron-proton 
interaction. This happens, in accordance with the sum rule (13), because the neutron and proton, which are close to 
each other in the phase-space, tend to exist in a bound not in a scattering state. And the n-p pairs, which form a 
deuteron, deplete the sample of n-p pairs and produce a dip of the correlation function at small relative momenta. 

The deuteron formation rate, which enters the sum rule (13), is computed with the deuteron wave function in the 
Hulthcn form 

, /<*/?(<* + /?)\V 2 e~ ar -e-P r 
Mr) - \2^W) r ' 

with a = 0.23 fm -1 and — 1.61 fm -1 [13]. Substituting the wave function (17) and the source function (4) into 
Eq. (9), we get 

A ■ - K a ^ a \S [*(2oro) - 2K((a + f3)r ) + K(2f3r )] , (18) 



rl («-/3) 2 



where 



K{x) 



2 2 f°° 

erfc(x) , erfc(x) = —= \ dte 

Jx 



In Fig. 3 we present the deuteron formation rate (18) as function of the source size parameter ro- As seen, Ad 
monotonously decreases when the source grows. 

In Figs. 4 and 5 we display the function ^(gmax), defined by Eq. (11), found for the singlet and triplet correlation 
functions presented in Figs. 1 and 2, respectively. Although S , (g max ) saturates at large (7 max for both the singlet 
and for the triplet correlation function, neither the sum rule (12) nor (13) is satisfied. The singlet S(q max ) does not 
vanish at large q max , as it should according to Eq. (12), while the triplet S(q niax ) is not negative, as required by 
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Eq. (13). However, comparing the numerical values of S(q m i iX ) to the corresponding deuteron formation rate, which 
sets the characteristic scale, one sees that with growing ro, the sum rule is violated less and less dramatically. It is 
not surprising as the asymptotic form of the wave function (14) and the s— wave approximation become then more 
accurate. The formula (11) shows that due to the factor q 2 even small deviations of the correlation function from 
unity at large q generate sizable contribution to S. Therefore, it can be expected that higher partial waves have to 
be taken into account to satisfy the sum rule. 



VI. HARD SPHERES 



We consider here a pair of the so-called hard spheres interacting via the potential which is zero at distances larger 
than a and it is infinite otherwise; a is the sphere's diameter. As well known, sec e.g. [14], scattering of the spheres 
is analytically tractable and the l~ th phase shift S\ is given by the equation 



tg<5. 



ni(qa) 

where ji(z) and n ; (z) is the spherical Bessel and Neumann function, respectively. The scattering amplitude equals 



(19) 



/( e) = I f> + 1)^0*6)1^, 



(20) 



where is the scattering angle and Pi(z) is the I— th Legendre polynomial. Substituting the amplitude (20) into the 
equation analogous to Eq. (14), we get the (asymptotic) wave function. Using this function and the Gaussian source 
(4), we get the correlation function of non-identical hard spheres as 



R(q) = l + —L r J2( 2l + l ) 



3i(qa) 



n, 



n 2 (qa)+jf(qa) 
(qa) cos(y) +ji(qa) sm(y 



(21) 



+ ( ni(qa) sin(y ) - j t (qa) cos 



drr e ir o cos(qr) ji(qr) 
drr e ir o sin(qr) ji(qr) 



3i(Qa) 



n 2 (l a )+3i(qa) ' 



For identical hard spheres, the wave function should be (anti-)symmetrized as ^ (0 q (r)±<^ q (— r)) , and the correlation 
function equals 



R(q) = 1 ± e- 4 ^ 2 + -pL- £(2Z + 1) 



ji(qa) 



n 2 (qa) +jf(qa) 



[1±(-1)' 



(22) 



(qa) cos (y) - ji(qa) sin( 



drr e 4r o cos (qr) ji(qr) 



-(ni(qa) sin (y) + 3i(q a ) cos (y)) J drr e ^ sin(qr) ji(qr) 



2r 2 q 2 y 



n 2 (qa)+jf(qa) 



[1± ("I)']- 



We note that there is a specific asymmetry of the signs of the corresponding terms in Eq. (21) and in Eq. (22). This 
happens because ((-1)' ± l)/2 = =p( — 1)'. 
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In Figs. 6-9 we show the correlation functions for non-identical spheres (21) and identical bosonic ones (22). The 
functions are computed numerically for the sphere diameter a = 1 fm and two values of the source size parameter 
tq = 4 fm and ro — 6 fm. We have considered rather large sources as the wave function, which is used to compute 
the correlation function, is of the asymptotic form (14). This requires r a. The functions S(q niax ) defined by 
Eq. (11), which correspond to the correlation functions shown in Figs. 6-9, are presented in Figs. 10-13. There are in 
Figs. 6-13 three families of curves representing: the s— wave approximation, the sum of partial waves with 1 = 0, 1, 2, 
and the sum of I = 0, 1, 2, 3, 4. As seen, going beyond the s— wave approximation minimally modifies the correlation 
functions shown in Figs. 6-9. It is rather expected as the correlation functions are significantly different from unity 
only at small momenta where contributions of higher partial waves are strongly suppressed. The higher partial waves 
are more important for the functions 5(g max ) but the integrals are still far from being saturated. According to the 
sum rule (8), the function S(q max ) should tend to zero for non-identical spheres and to n 3 T> r (0) = (y / 7r/2ro) 3 for 
identical ones when <7 max — > oo. Such a behavior is not observed in our calculations and taking into account the 
partial waves of higher I does not improve the situation. We see two possible explanations of the problem. 1) The sum 
rule is in principle fulfilled but one should go to much larger momenta q max to saturate the integral (11). However, 
higher relative momenta require taking into account more and more partial waves. If this is indeed the case, the sum 
rule is formally correct but useless because Eq. (1), which is a starting point of the sum rule derivation, assumes that 
the relative momentum of particles q is much smaller than the typical particle momentum. 2) The integral (11) is 
divergent as g max goes to infinity. Then, the sum rule is simply meaningless. We discuss the second possibility in 
more detail in the next section where the Coulomb interaction is studied. 



VII. COULOMB INTERACTION 



As well known, the Coulomb problem is exactly solvable within the non-relativistic quantum mechanics [14]. The 
exact wave function of two non-identical particles interacting due to the repulsive Coulomb force is given as 

q (r) = T(l + i^) e^ z ' 2 F(-i^, 1, iq V ) , (23) 

where A = /xe 2 /87r with /x being the reduced mass of the two particles and ±e is the charge of each of them; F 
denotes the hypergeomctric confluent function, and r\ is the parabolic coordinate (see below). The wave function 
for the attractive interaction is obtained from (23) by means of the substitution A — ► —A. When one deals with 
identical particles, the wave function ^> q (r) should be replaced by its (anti-)symmetrized form. The modulus of the 
wave function (23) equals 

|0 q (r)| 2 = G(q) \F(-i±l,iqr,)\ 2 , 

where G(q) is the so-called Gamov factor defined as 

_,, „ 27rA 1 ,„ 

G(q) = . 24 

q cxp( 2 ^)-l 

As seen, the modulus of the wave function of non-identical particles solely depends on the parabolic coordinate rj. 
Therefore, it is natural to calculate the Coulomb correlation function in the parabolic coordinates: r\ = r — z, £ = r + z 
and 4> which is the azimuthal angle. Then, the correlation function computed for the Gaussian source function (4) 
equals 

m = £jk[ d ^<-^) l^-VMrf, (25) 

where the integration over £ has been performed. It should be noted here that, in contrast to the neutron-proton 
and hard sphere cases, the Coulomb correlation function is calculated with the exact wave function not with the 
asymptotic form of it. 

The modulus of the (anti-)symmetrized Coulomb wave function equals 

I^WI 2 = \ G{q) \\F(-i-, 1, iq V )\ 2 + \F(-i-, 1, iq£)\ 2 ± 2Re(e^ F(-i-, 1, iqrf) F*(-i-, 1, itf)) 
z l q q v 9 Q ' 

and the correlation function analogous to (25) is 
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R(q) = 




(26) 



In Figs. 14 and 15 we demonstrate the correlation function of non-identical particles and of identical pions, respec- 
tively, which are given by Eqs. (25,26). In the case of non-identical particles, we have taken the reduced mass, which 
enters the parameter A, to be equal to that of identical pions. Therefore, the two systems differ only due to the effects 
of quantum statistics. In both cases the Coulomb interaction is repulsive. While the correlation of non-identical 
particles is everywhere negative, the correlation of identical pions is negative at small q due to the Coulomb repulsion 
but it is positive for larger q due to the effect of quantum statistics. 

Figs. 16, 17, and 18 show the functions 5(<7 max ) defined by Eq. (11) which are computed for the correlation functions 
presented in Figs. 14 and 15, respectively. According to the sum rule (8), S(q m:ix ) should vanish for sufficiently large 
<Zmax in the case of non-identical particles which repel each other, and 5(g m ax) should tend to 7r 3 T> r (0) for identical 
pions. As seen in Figs. 16, 17, and 18, the functions S(q ma ^) do not seem to saturate at large <7 max and the sum rule 
is badly violated. 

What is wrong here? The derivation of the sum rule (8) implicitly assumes that the integral in the l.h.s. exists i.e. it 
is convergent. Otherwise interchanging of the integrations over q and over r, which leads to Eq. (5), is mathematically 
illegal. Unfortunately, the Coulomb correlation functions appear to decay too slowly with q, and consequently the 
integral in the left hand side of Eq. (8) diverges. 

To clarify how the integral diverges, one has to find the asymptotics of the correlation function at large q. This 
is a difficult problem as one has to determine a behavior of the wave function (23) at large q for any r, and then to 
perform the integration over r. Up to our knowledge, the problem of large q asymptotics of the Coulomb correlation 
function has not been satisfactory solved although it has been discussed in several papers [15-19]. We have not found 
a complete solution of the problem but our rather tedious analysis, which uses the analytic approximate expressions 
of the hypergeometric confluent function at small and large distances, suggests, in agreement with [15-19], that 
(R(q) — 1) ~ l/<7 2 when q — > oo. Then, the integral in the left hand side of Eq. (8) linearly diverges, and the sum 
rule does not hold. Wc note that the Gamov factor (24), which represents a zero size source and decays as 1/q at 
large q, leads to the quadratic divergence of the integral (8). We also note here that the asymptotics 1/q 2 of the 
correlation function does not have much to do with the well known classical limit of the correlation function [15-19]. 
Since the large q limit of the correlation function corresponds to the small separation of the charged particles, which 
at sufficiently large q is smaller than the de Broglie wave length, the classical approximation breaks down. 



The sum rule (8) provides a rigorous constraint on the correlation function if the momentum integral in Eq. (8) 
exists. The rule is trivially satisfied by the correlation function of non-interaction particles. The model calculations 
of the n-p correlation function in the s— wave approximation fail to fulfill the sum rule. One suspects that the 
approximation, which is sufficient to properly describe a general shape of the correlation function, distorts its tail. 
Since even small deviations of the correlation function from unity at large q generate a sizable contribution to the 
sum-rule integral, it seems reasonable to expect that the higher partial waves have to be included to comply with 
the sum rule. Paradoxically, when the higher partial waves are taken into account for the interacting hard spheres 
the situation is not improved but it gets worse. This suggests that either one should go to really large momenta to 
saturate the sum-rule integral or the integral in Eq. (8) is divergent. In the case of Coulomb repelling interaction, 
we certainly deal with the second option. Due to the strong electrostatic repulsion at small distances, the correlation 
function decays too slowly at large momenta, and consequently the sum- rule integral does not exist. 

Being rather useless, the sum rule explains some qualitative features of the correlation function. In particular, it 
shows that, in spite of the attractive interaction, the correlation can be negative as it happens in the triplet state of 
neutron-proton pair and in the pA system. 



We are very grateful to Wojcicch Florkowski for a generous help in performing numerical calculations, and to David 
Brown, Cristina Manuel, Richard Lcdnicky, and Scott Pratt for fruitful discussions and correspondence. 
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FIG. 1. The singlet correlation function of neutron and proton computed for three values of the source size parameter ro- 
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FIG. 2. The triplet correlation function of neutron and proton computed for three values of the source size parameter ro. 




2 4 6 8 10 



FIG. 3. The deuteron formation rate as a function of the source size parameter ro. 
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FIG. 4. The function <S(g ma x) for to the singlet correlation function of neutron and proton for three values of the source size 
parameter ro- To set the scale the corresponding values of the deuteron formation rate are given. 
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FIG. 5. The function S , (g m ax) for to the triplet correlation function of neutron and proton for three values of the source size 
parameter ro- 
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FIG. 6. The correlation function of non-identical hard spheres for the source size parameter ro = 4 fm. 
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FIG. 7. The correlation function of non-identical hard spheres for the source size parameter ro = 6 fm. 
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FIG. 8. The correlation function of identical hard spheres for the source size parameter ro = 4 fm. 
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FIG. 9. The correlation function of identical hard spheres for the source size parameter ro = 6 fm. 
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FIG. 11. The function S(q ma ^) corresponding to the correlation function of non-identical hard spheres for the source size 
parameter r — 6 fm. 
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FIG. 12. The function S^gmax) corresponding to the correlation function of identical hard spheres for the source size parameter 
ro = 4 fm. 




FIG. 13. The function S , (g m ax) corresponding to the correlation function of identical hard spheres for the source size parameter 
ro = 6 fm. 
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FIG. 14. The Coulomb correlation function of non-identical repelling particles for three values of the source size parameter 
ro- 
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FIG. 15. The Coulomb correlation function of identical pions for three values of the source size parameter ro. 
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FIG. 16. The function S'(q max ) corresponding to the Coulomb correlation function of non-identical repelling particles for 
three values of the source size parameter r . 




FIG. 17. The function S , ((j , ma x) corresponding to the Coulomb correlation function of identical pions for the source size 
parameter ro = 1 fm. The dotted line represents the value of 7r 3 V r (0). 
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FIG. 18. The function S(q mB ^) corresponding to the Coulomb correlation function of identical pions for the source size 
parameter ro equal 3 fm and 5 fm. The dotted lines represents the respective values of ir 3 T> r (0). 
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